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ABSTRACT 

We present measurements of the angular correlation function of galaxies selected from 
a Bj ^ 23.5 multicolour survey of two 5° x 5° fields located at high galactic latitudes. 
The galaxy catalogue of ~ 4 x 10^ galaxies is comparable in size to catalogues used to 
determine the galaxy correlation function at low-redshift. Measurements of the z ~ 0.4 
correlation function at large angular scales show no evidence for a break from a power 
law though our results are not inconsistent with a break at ^ 15/i~^Mpc. Despite 
the large fields-of-view, there are large discrepancies between the measurements of 
the correlation function in each field, possibly due to dwarf galaxies within z ^ 0.11 
clusters near the South Galactic Pole. 

Colour selection is used to study the clustering of galaxies z to z 0.4. The 
galaxy correlation function is found to strongly depend on colour with red galaxies 
more strongly clustered than blue galaxies by a factor of ^ 5 at small scales. The slope 
of the correlation function is also found to vary with colour with 7 '-^ 1.8 for red galaxies 
while 7 ~ 1.5 for blue galaxies. The clustering of red galaxies is consistently strong 
over the entire magnitude range studied though there are large variations between 
the two fields. The clustering of blue galaxies is extremely weak over the observed 
magnitude range with clustering consistent with vq ~ 2/i~^Mpc. This is weaker than 
the clustering of late- type galaxies in the local Universe and suggests galaxy clustering 
is more strongly correlated with colour than morphology. This may also be the first 
detection of a substantial low redshift galaxy population with clustering properties 
similar to faint blue galaxies. 

Key words: (cosmology:) large-scale structure of Universe - galaxies: evolution. 



1 INTRODUCTION 

The galaxy two-point correlation function is commonly used 
to measure the structure of the galaxy environment from 
high redshift until the present epoch. The clustering prop- 
erties of galaxies in the local Universe are well measured by 
large representative surveys of the galaxy population (Mad- 
dox, Efstathiou & Sutherland 1996). Catalogues of galaxies 
selected by morphology show large variations of the galaxy 
correlation function with late type galaxies having consid- 
erably weaker clustering than early type galaxies (Davis & 
Geller 1976, Loveday et al. 1995). 

The results from studies of galaxies with fainter ap- 
parent magnitudes and higher redshifts are less conclusive. 
Pencil-beam surveys with CCDs and photographic plates 
from 4m telescopes have measured the amplitude of the 
B > 22 correlation function; however, estimates vary by 
Z 100% (Infante & Pritchet 1995). Also, while B > 22 sur- 
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veys show evidence for a rapid decline of the amplitude of 
the correlation function (Efstathiou et al. 1991, Infante & 
Pritchet 1995, Roche et al. 1996), I band imaging surveys 
to similar depths show no evidence for a rapid decrease of 
the correlation function amplitude (Postman et al. 1998). 

The small areas of previous studies of the faint galaxy 
correlation function are a possible source of the discrepancy. 
Large individual structures and voids in the Universe could 
bias estimates of the correlation function if the field-of-view 
of the survey is small. The use of single band data to se- 
lect catalogues of galaxies could suffer from biases as the 
morphological mix of galaxies will change as a function of 
limiting magnitude. It is probable that the differing ampli- 
tudes of the B and I band correlation functions are due to 
faint B band data being dominated by weakly clustered blue 
galaxies (Efstathiou et al. 1991) while the / band data has 
a larger fraction of early type galaxies. 

In this paper, we use a Bj ~ 23.5 multicolour catalogue 
of galaxies derived from two 5° x 5° fields to measure the 
clustering properties of faint galaxies. Section 2 discusses the 
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observations and data reduction used to produce the galaxy 
catalogue. The method used to determine the angular and 
spatial correlation functions is described in Sections 3 and 4. 
Estimates of the correlation function at large angular scales 
are presented in Section 5. Section 6 discusses the angular 
correlation function as a function of limiting magnitude for 
the U, Bj, Rf and I bands. We show that it is impossible 
to use a single model to describe the observed clustering of 
galaxies as different populations of galaxies are measured by 
each band as a function of limiting magnitude. In Section 
7, we use colour selection to measure galaxies with similar 
stellar populations over a range of redshifts. Our estimates 
of the spatial correlation function indicate clustering is more 
strongly correlated with colour than morphology. Our main 
conclusions are summarised in Section 8. 



2 THE GALAXY CATALOGUE 

The image data consists of U, Bj, Rf and / band digi- 
tally coadded SuperCOSMOS scans of photographic plates 
of the South Galactic Pole Field (SGP) and UK Schmidt 
field 855 (F855). The field centres, number of photographic 
plates and limiting magnitudes for each field and band are 
listed in Table 0. The individual photographic plates were 
obtained with the 1.25m UK Schmidt Telescope between 
1977 and 1998 for various programmes. The field-of-view is 
6.4° X 6.4° though vignetting is greater than 0.1 magnitudes 
at angular scales more than ~ 3.35° from the optical axis 
(Tritton 1983). 

The individual plates were digitised with the Super- 
COSMOS facility at Royal Observatory Edinburgh (MiUer 
et al. 1991). The 10/im pixels and 15/im resolution of Super- 
COSMOS corresponds to a pixel scale of 0.67" and allows 
good sampling of images obtained in typical seeing (~ 2.5"). 
The plate scans were converted from transmission to ap- 
proximate intensity space before stacking. All plates were 
sky-limited so faint objects are in the linear regime of the 
plate response. 

Before coadding the individual plate scans, the back- 
ground was subtracted with a 160 x 160 pixel median filter. 
Without background subtraction, robust bad pixel rejection 
is not possible resulting in contamination of the catalogue. 
To reduce the computation time required to subtract the 
background, the background was calculated for 32 x 32 pixel 
regions rather than for each pixel and only every 4th pixel 
in the 160 x 160 filter box was used to determine the me- 
dian. Before coadding the scans, the noise distribution of the 
individual plates was measured and found to be fitted ac- 
curately by a Gaussian curve. If the noise distribution was 
significantly skewed, large systematic errors would propa- 
gate into the coadded plate scans. 

As faint objects are in the linear regime of the plate 
response and the noise distribution is approximately Gaus- 
sian, the IRAF0 task imcombine can be effectively used to 
coadd the data. An average of the plates weighted by the 

t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



background noise and min-max bad pixel rejection is used 
to provide high signal-to-noise and robust bad pixel rejec- 
tion. Without robust bad pixel rejection, there is significant 
contamination of the catalogue by plate-flaws, satellite trails 
and asteroids. For bands with 5 or less plates, the number of 
pixels rejected by the bad pixel rejection algorithm results 
in the coaddition effectively being a median. For the SGP U 
band data, where only 2 of the 3 plates are of high quality, 
a weighted average with no bad pixel rejection is used and 
objects are correlated with the B,j, Rf and / catalogues 
to remove contamination. The U band data is complete to 
[/ = 22 for objects with B,j detections but incompleteness 
does effect the catalogue at (7 > 21.5 for faint galaxies with 
U-Bj < -1.5. 

Object detection, instrumental photometry and faint 
object star-galaxy classifications are determined using SEx- 
tractor (Berlin & Arnouts 1996). In the non- linear regime 
of the plate response where SExtractor classifications are 
unreliable, objects are classified as stellar if they are within 
the stellar locus in peak intensity versus magnitude space 
and area versus magnitude space. The reliability of the ob- 
ject classifications is improved by using object classifications 
in multiple bands where possible. Objects with SExtrac- 
tor classification scores greater than 0.7 in 3 bands, 0.75 
in 2 bands or 0.85 in a single band are classified as stel- 
lar objects. Comparison of the object classification with the 
Bj < 22 and i? < 24 spectroscopic samples of CoUess et al. 
(1990) and Glazebrook et al. (1995) show the star-galaxy 
separation to be correct for ~ 85% of Bj < 23 objects and 
contamination of the galaxy catalogue is ^ 5%. At Bj < 20, 
comparison of object classifications with CCD images, eye- 
ball classifications and colour-colour diagrams indicates the 
reliability of the star-galaxy separation is <: 95%. 

Photometric calibration of the data was obtained with 
B, V , R and I band CCD data from the Siding Spring 40- 
inch provided by Bruce Peterson, B and R band CCD imag- 
ing with the Anglo- Australian Telescope provided by Bryn 
Jones and U band photometry from Groom et al. (1999) and 
Osmer et al. (1998). Galaxy and faint stellar photometry is 
obtained by adding a zero-point value to the instrumental 
magnitudes while bright stellar photometry is determined by 
a polynomial fit between the instrumental magnitudes and 
calibrated CCD magnitudes. All magnitudes are corrected 
for vignetting using a fit to the vignetting as a function 
of radius plot from Tritton (1983). Comparison of the pho- 
tometry with Boyle, Shanks & Croom (1995), Caldwell & 
Schechter (1996), Croom et al. (1999), Osmer et al. (1998) 
and Patch B of the ESO Imaging Survey (Prandoni et al. 
1999) limits zero point errors to less than 0.05 magnitudes. 
Comparison of median colours as a function of limiting mag- 
nitude and galaxy loci in colour-colour diagrams show no 
significant offsets between the two fields. 

Magnitude estimates for galaxies are corrected for dust 
extinction using the dust maps of Schlegel, Finkbeiner & 
Davis (1998). Figure |l| shows that there are significant vari- 
ations in the dust extinction in F855 which could introduce 
spurious large-scale structure into the catalogue. In contrast, 
the dust extinction in the SGP is restricted to the range 
0.01 < E{B -V) < 0.03 which is a comparable to the 0.028 
magnitude error estimate of the dust maps. We therefore 
use the dust map estimates of the extinction for F855 while 



© 1999 RAS, MNRAS 000, 



The Clustering of Colour Selected Galaxies 3 



Table 1. Image data parameters 
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Figure 1. A plot of the dust extinction across the SGP and F855 
fields using E{B — V) estimates from Schlegel, Finkbeiner & Davis 
(1998). The dust extinction in F855 {I ~ 45°) is significantly 
larger and shows more structure than the SGP. The extinction in 
F855, if left uncorrected, would introduce spurious structure on 
large angular scales. 



using a constant value of E{B — V) = 0.015 to correct for 
dust extinction in the SGP. 

Completeness limits for the catalogue are determined 
with comparison to deeper CCD images (where available) 
and number counts as a function of limiting magnitude 
(Figure |^. The number counts for the two fields are in 
good agreement with each other in all 4 bands. However, 
at Bj < 19 and Bj > 21 the number counts are slightly 
higher than those measured by Maddox et al. (1990) and 
Infante & Pritchet (1992). The resulting completeness lim- 
its for the four bands are U — 21.5, Bj = 23.5, -Rf = 22 and 
7 = 20 for the SGP and U = 22, Bj = 23.5, Rf = 22.5 and 
/ = 20 for F855. Overestimates of the completeness limits 
would result in biases in estimates of the correlation func- 
tion as vignetting would introduce spurious structure into 
the catalogue. 

Contamination of the catalogue by spurious objects bi- 
ases estimates of the galaxy correlation function and these 
have to be removed. Most plate-flaws, asteroids and satellite 
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Figure 2. A plot of the galaxy number counts as a function of 
magnitude (after correction for dust extinction). Data from both 
fields is compared with number counts from the APM survey 
(Maddox et al. 1990) and Infante & Pritchet (1992). / band counts 
have been moved 1 magnitude to the left and U band counts have 
been moved 2 magnitudes to the right for clarity. 



trails are effectively removed by the bad-pixel rejection al- 
gorithm. The plate edges, the step wedges, globular clusters, 
bright galaxies and the internal reflections of y < 5 stars 
are manually removed from the catalogue. The majority of 
the remaining contaminating objects are difi'raction spikes 
and halos around bright stars and these are removed with 
a program that drills regions around bright stars. Spurious 
structure may also be introduced into the catalogue at an- 
gular scales more than 3.35° from the plate centre where 
number counts are altered by > 10% during the vignetting 
correction. To prevent this, the corners of the fields have 
been excluded from the catalogue. The resultiiig geometry 
of the 5° X 5° field-of-view is shown in Figure H. The final 
catalogues of galaxies contain ~ 2 x 10° galaxies per field in 
both Bj and Rf- 



3 ESTIMATION OF THE ANGULAR 
CORRELATION FUNCTION 

The galaxy two-point correlation function, a; (61), measures 
the mean excess surface density of pairs at angular separa- 
tion 6 compared with the expected number of pairs if galax- 
ies were randomly distributed. The most commonly used 
estimator of the angular correlation function is 

DD 
DR 

where DD and DR are the number of galaxy-galaxy and 
galaxy-random object pairs at angular separations 9 ± 59. 
The random objects are typically copies of real objects dis- 
tributed randomly across the field-of-view. To reduce errors, 
multiple random copies of each object can be made and the 
estimate of DR renormalised. However, the DD/DR esti- 



cbi9) 



(1) 
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Figure 3. A flux weighted map of Rp < 22 galaxies in tlie SGP. 
The dark blocks are regions which have been removed from the 
catalogue. The corners of both fields have been removed to pre- 
vent objects with significant vignetting entering the catalogue. 
Some of the Abell clusters in the field can be seen as overdcnse 
regions. 



mator is subject to first-order errors in the galaxy density 
contrast (Hamilton 1993) making it unsuitable for measur- 
ing weak clustering at large angular scales. We therefore use 
an estimator with lower variance than DD/DR, the estima- 
tor 

DD - 2DR + RR 



RR 



(2) 



(Landy & Szalay 1993) where RR is the number of random- 
random object pairs at angular separations 6±56. The value 
of u}{9) is determined for different angular scale bins with 
the values of DD, DR and RR being determined with pairs 
of individual objects at small angular scales and weighted 
pairs of cells containing multiple objects at large angular 
scales. The use of cells to determine the angular correlation 
function reduces the computational time required to several 
hours and does not introduce significant errors. 

The estimator of the angular correlation function satis- 
fies the integral constraint, 



(3) 



(Groth and Peebles 1977), resulting in an underestimate of 
the angular correlation function. To remove this bias from 
the correlation function, the term 



(4) 



is added to the estimate of the correlation function. The 
term, Lu{9)n does require an assumption of the form of the 
correlation function to correctly estimate the value of corre- 
lation function. However, previous work with smaller fields 
of view and at brighter limiting magnitudes shows that the 



angular correlation is well approximated by a power law at 
angular scales less than 1° (Maddox, Efstathiou & Suther- 
land 1996). 

A further source of bias in estimates of the correlation 
function is contamination of the catalogue by randomly dis- 
tributed objects such as stars. If the fraction of the cata- 
logue contaminated by stars is /, the estimate of the cor- 
relation function is reduced by a factor of (1 — /)^ at all 
angular separations. As mentioned previously, comparison 
with spectroscopic samples indicates the star galaxy separa- 
tions if ^ 85% reliable at magnitude brighter than Bj < 22. 
At fainter magnitudes, the contamination of the catalogue 
by stars is negligible as galaxies consist more than ^ 80% 
of the object number counts at Bj > 23 (Glazebrook et al. 
1995). 



4 MODELLING THE SPATIAL 
CORRELATION FUNCTION 

For this work, we assume that the spatial correlation func- 
tion is a power law of the form. 



+ 



(5) 



(Efstathiou et al. 1991) where r is the spatial separation in 
physical coordinates, z is the redshift and ro, 7 and e are con- 
stants. If e = the clustering is fixed in physical coordinates 
while if e = 7 — 3 the clustering is fixed in comoving coordi- 
nates. This parameterisation of the evolution of the spatial 
correlation function is not valid at all redshifts but is a good 
approximation at 2 < 1 (Baugh et al. 1999). For galaxies 
selected with images in a single broadband, typical values 
of the parameters of ^(''i-^) s-r^ ''0 ~ 5/i~^Mpc, 7 ~ 1.7 

(Maddox, Efstathiou & Sutherland 1996) and 1. 

For a power law spatial correlation function, the result- 
ing angular correlation function is a power law with 



,r[(7-i)/2] 



,7fl(l — >) 



(6) 



r(7/2) 

(Baugh & Efstathiou 1993) where B is a constant. The value 
of B is given by 



B 



where 



dN{z) 
dz 



dz 



dN{z] 
dz 



dz 



aiz) 



dz 
dx 



1-7 



Fix)il + z)-'-^+'-'\ 



(7) 



(8) 



X is the coordinate distance at redshift z, dN/dz is the num- 
ber of galaxies per unit redshift detected by the survey and 



F{xf ^l + nR{Hox/cf. 
The value of x is given by 

_ c r 1 

Ho J, E{z) 
where 

E{z) = v/nA/(l + 2)^ + !^j?(l + 2)2 + f^A. 



(9) 



(10) 



(11) 



For values of z < 0.5, the value of B is more strongly de- 
pendent on ro, 7 and the galaxy redshift distribution than 
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the cosmological model. This is not unexpected as the value 
of X for z — 0.4 varies by less than 10% between Qm = 1 
and Qm = 0.2 models of the Universe. We therefore only 
use a single cosmological model with Ho — 75kms~^Mpc~^, 
= 0.2 and = 0. 

There are two approaches to modelling the galaxy red- 
shift distribution which are often used in the literature. The 
first approach uses an accurate description of the local lu- 
minosity function of different galaxy types, plus models of 
galaxy evolution and fc-corrections for each type, and at- 
tempts to model the observed number counts and redshift 
distribution where data is available (i.e. Roche et al. 1996). 
This model includes the physics of galaxy evolution, how- 
ever, it contains large numbers of free parameters and dif- 
ferent models can readily reproduce the observed number 
counts. The second approach, which is applied to our single 
band imaging data, assumes a functional form for the red- 
shift distribution and uses the observed number counts and 
redshift surveys to constrain the model (Baugh & Efstathiou 
1993). This approach produces a good model of the redshift 
distribution but contains no physics of galaxy evolution and 
is limited by the depth of redshift surveys. 

To model the redshift distribution, we use the model 
of Baugh & Efstathiou (1993). For a complete sample of 
galaxies brighter than magnitude m, the number of galaxies 
detected per unit redshift is given by 



Survey Depth 



dN(z) ^, ,2 

= A(m)z exp 



dz 



where 



A{m) 



1.4122„(m) 



3/2N 



1.062N{m)n 

Zmim) 



(12) 



(13) 



N{m) is the galaxy number counts, Q is the survey area and 
Zm is the median redshift of the galaxy sample. As the value 
of B depends on the galaxy redshift distribution and not the 
galaxy number counts, the most important free parameter 
is Zrri which is a function of band and limiting magnitude. 

The median redshift as a function of limiting magni- 
tude is obtained from a polynomial fit to median redshifts 
derived from galaxy redshift surveys. At low redshift this 
is derived from the local galaxy luminosity function while 
at higher redshifts the median redshifts determined directly 
from the survey data are used. Local luminosity functions 
are assumed to be Schechter functions with a — —1.00, 



Mb J =-19.5-51og/i 
Loveday et al. (1992) and 



-20.5 - 51og/i. 



(14) 



(15) 



Changing the value of a to —1.3 reduces the estimate of the 
median redshift by ~ 15% resulting in a similar decrease of 
the estimate of ro. Galaxy fc-corrections are approximated 
by ksiz) — 2z and kR{z) = 0.5z. At z ^ 0.2, median red- 
shifts as a function of survey depth are derived from the 
redshift surveys of CoUess et al. (1990), Glazebrook et al. 
(1995), Lin et al. (1999) and Munn et al. (1997). / band me- 
dian redshifts are the same as Postman et al. (1998) which 
were derived from the CFHT redshift survey (Lilly et al. 
1995). To determine the U band median redshift, we use 
U — Bj ~ 0.3, the median colour of Bj < 21 galaxies in the 
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Figure 4. A plot of the model median redshift as a function of 
survey depth for U, Bj, Rp and / bands. The median redshift for 
the U band is derived using a median galaxy colour of U — Bj ~ 
0.3. Magnitude limits for the SGP and F855 are shown with the 
large filled and open circles respectively. 



F855 field. The functions for median redshift as a function 
of limiting magnitude are shown in Figure ^. 



5 THE ANGULAR CORRELATION 

FUNCTION AT LARGE ANGULAR SCALES 

The 5° X 5° field-of-view of each field allows the measure- 
ment of the faint galaxy correlation function at large angular 
scales. At Bj ~ 23.5 and Rf ~ 22, the median redshift of 
the data is 2; ~ 0.4 and a break in the spatial correlation 
function at ~ 15/i~^Mpc in comoving coordinates (Maddox, 
Efstathiou & Sutherland 1996) corresponds to an angular 
scale of ~ 1°. 

It is possible to measure the angular correlation for each 
field to ~ 5° but at large angular scales the estimate of the 
correlation function will be dominated by individual struc- 
tures. To determine the range of angular scales where the 
correlation function is representative, we compare the esti- 
mates of the angular correlation function for the SGP and 
F855 fields. This provides a more reliable estimate than sub- 
samples of the data as each subsample would have smaller 
field-of-view than the original data, resulting in overesti- 
mates of the errors at large angular scales. 

The Bj < 23.5 and Rf < 22 correlation functions for 
each field are shown in Figures ^ and |^. As the integral con- 
straint depends of 7, the data has been fitted with power 
laws with fixed 7 to allow comparison of the 2 fields es- 
timates of uj{0). For both fields the integral constraint is 
less than the amplitude of the correlation function at 1°. A 
power law matches the data well and there is no evidence 
of a break from a power law on all angular scales. However, 
on angular scales ^ 2°, the estimates of u}{9) for both fields 
are within 2(j of 0. It is therefore possible that the break in 
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17.0<Bj<23.5 Galaxy Angular Correlation Function 
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Figure 5. The Bj < 23.5 angular correlation function for the 
SGP and F855 fields. Power law fits to the SGP and F855 data 
are shown with solid and dashed lines respectively. The data in 
each field is fitted well by a power law and most data points are 
within 2cr of each other. There is no evidence of a break in the 
correlation function on any of the scales measured though lo{9) at 
> 1° is only ~ 2a more than 0. It is also possible both correlation 
functions could be biased by large structures on scales > 1°. 



15.5<Rp<22.0 Galaxy Angular Correlation Function 




8 (Degrees) 

Figure 6. The Rp < 22 angular correlation function for the 
SGP and F855 fields. Power law fits to the SGP and F855 data 
are shown with solid and dashed lines respectively. The data for 
each field is fitted well by a power law and there is no evidence 
of a break in the correlation function on scales less than ~ 2°. 
However, unlike the Bj < 23.5 data, there is an offset between 
the data in the two fields. 

the correlation function is present at 15/i~^Mpc in comoving 
coordinates but can not be detected witli this dataset. 



6 THE CORRELATION FUNCTION AS A 
FUNCTION OF LIMITING MAGNITUDE 

The large fields-of-view used for this study reduce the errors 
associated with large structures along the line-of-sight. Also, 
the input catalogue of up to 2.5 x 10^ galaxies results in small 
random errors in estimates of u>{9). Previous measures of the 
correlation function in U, Bj and Rf have been restricted to 
fields-of-view less than ~ 2° x 2° (Infante & Pritchet 1995). 



At the present time, only the / band survey by Postman et 
al. (1998) has a comparable field-of-view with greater depth 
than this work. 

Previous measurements of 7 indicate that it may vary 
with band and survey depth (Infante & Pritchet 1995, Post- 
man et al. 1998) though most values in the literature are 
between 1.6 and 1.8. The consistent data reduction method 
used for this work should allow the accurate comparison of 
7 as a function of band and survey depth. 

Estimates of the value of a;(l') and Tderived from power 
law fits to the data are shown in Table H. The amplitude of 
the correlation function is determined at 1' rather than 1° 
as the amplitude and 7 are not independent and estimates 
of w(l°) depend strongly on 7. The values of 7 show a weak 
trend towards smaller values with magnitude and bluer sur- 
vey bands. As there is a correlation between galaxy colour 
and morphology (see Figure |l^ and late type galaxies have 
shallower values of 7 than early type galaxies, this trend is 
not unexpected. 

While the value of 7 does not differ by more then ~ 2a 
between the two fields, the amplitude of the correlation func- 
tion varies by ~ 100% at bright magnitudes. While varia- 
tions of amplitude could be caused by zero point errors, the 
error required is approximately 0.4 magnitudes in all bands 
at bright magnitudes with it decreasing to ~ at Bj ~ 23. 
This is inconsistent with the photometric calibration, the 
galaxy number counts in Figure ^ and the colour-colour di- 
agrams. A systematic error could be present in the data but 
it seems unlikely that it could effect lo{1') without causing 
large variations of 7. 

If the values of do not have significant system- 

atic errors, a possible cause of the variations is that the 
two fields measure different populations of galaxies with dif- 
ferent clustering properties. The SGP, which has stronger 
clustering than F855, contains "sheets" perpendicular to 
the line-of-sight (Broadhurst et al. 1990) and smaller struc- 
tures including 22 Abell clusters (Abell, Corwin & Olowin 
1989). Of these clusters, 11 appear to be associated with a 
structure at z ~ 0.11 identified by Broadhurst et al. (1990). 
Mbj ~ —19.5 (M*) galaxies at z ~ 0.11 have an appar- 
ent magnitude of B,j ~ 18 and this population of galeixies 
is too small to significantly bias estimates of the Bj > 20 
correlation function. However, it is possible that the cluster 
population of Mb , ^ — 16 galaxies could bias estimates of 
the correlation function if they are a significant fraction of 
the observed galaxy number counts. 

To test if the z — 0.11 clusters do significantly effect 
the correlation function, the correlation function has been 
determined with 2° x 2° (~ 10/i"^Mpc x 10/i"^Mpc) re- 
gions surrounding the clusters removed from the data. As 
the clusters are not uniformly distributed across the field- 
of-view, the size of the catalogue is reduced by ~ 60%. Ta- 
ble ^ lists the amplitude of the Bj and Rf band correlation 
functions for the SGP field without the clusters and F855 
for comparison. The amplitude of the correlation function 
has decreased significantly compared with the original esti- 
mates for the SGP. While most noticeable at bright mag- 
nitudes, the effect is also significant at fainter magnitudes 
where the contribution from nearby clusters might be ex- 
pected to be small. This is consistent with Mbj ^ —16 
galaxies within the clusters significantly effecting estimates 
of the faint galaxy correlation function. However, details of 
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Table 2. Measured parameters for ii>(9) determined from power law fits to data between 10" and 0.3°. 

Field SGP F855 



Magnitude Range 
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Table 3. Measured parameters for ijj{0) determined from power law fits to data between 10" and 0.3°. Clusters at 2 ~ 0.11 have been 
removed from the SGP sample to reduce the effect dwarf galaxies may have on the estimate of the correlation function. 

Field SGP (no z = 0.11 clusters) F855 



Mag 


nitude Range 
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the relationship between clusters and the observed cluster- 
ing of faint galaxies will be explored in a later paper. 

For the following discussion of the amplitude as a func- 
tion of limiting magnitude, the estimates of the SGP correla- 
tion function including the z = 0.11 clusters are used. While 
it is probable relatively nearby clusters do effect estimates 
of the correlation function, there is no clear justification for 
excluding them. Also, excluding the regions surrounding the 
clusters significantly reduces the number of galaxy pairs used 
to determine the correlation function, significantly increas- 
ing random errors. 

Figures ^ to ^ plot the I, Rp, Bj and U band angu- 
lar correlation function amplitudes as a function of limiting 
magnitude. The amplitude of the correlation function has 
been determined with fixed values of 7 to reduce the depen- 



dence of the amplitude on 7. As the values of the angular 
correlation function differ significantly between each field, 
no attempt has been made to fit the data. Instead, a model 
has been plotted with clustering fixed in comoving coor- 
dinates and ro — S/i^'^Mpc. The value of ro is similar to 
values derived by Maddox, Efstathiou & Sutherland (1996) 
and Postman et al. (1998). 

Figure ^ shows the measured amplitude of the I band 
correlation function for the SGP, F855 and Postman et al. 
(1998). The measured clustering in the SGP is in good 
agreement with Postman et al. (1998), while F855 measures 
slightly weaker clustering. The model of the clustering is also 
a reasonable estimate of the observed clustering across the 
magnitude range observed. Figures ^ to |l^ show the Rp, 
Bj and U band clustering to be significantly weaker than 
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Figure 7. The / band correlation function amplitude. The mea- 
surements of the amplitude are consistent with the data of Post- 
man et al. (1998) which is shown with crosses. The data from 
Postman et al. (1998) has been corrected for the assumption that 
7 = 1.7. 
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Figure 8. The Rp band correlation function amplitude. While 
the redshift range is similar to the / band data, the galaxy clus- 
tering is significantly weaker. 



the / band clustering for galaxies with a similar range of 
redshifts. Comparison with the results of Infante & Pritchet 
(1995) show the SGP has similar clustering in Bj and weaker 
clustering in Rp- Infante & Pritchet (1995) may measure 
stronger clustering than F855 as their field is located in the 
NGP which has similar large scale structures as the SGP 
(Broadhurst et al. 1990). 

While the ro = 5.0/i~^Mpc model is a good fit to the 
SGP Rf < 20 and Bj < 21 data, at fainter magnitudes a 
rapid decline in the amplitude of the correlation function is 
observed. The rapid decline of the faint correlation function 
has been observed previously by Efstathiou et al. (1991), In- 
fante & Pritchet (1995) and Roche et al. (1996). Obviously, 
the Rf and Bj band data are sampling a different popula- 
tion of galaxies to the / band sample. To be consistent with 
galaxy redshift surveys, the Bj ~ 23 galaxies must be domi- 
nated by a population of weakly (ro ~ 2.0/i~^Mpc) clustered 
galaxies at « ~ 0.4 (Efstathiou et al. 1991, Efstathiou 1995). 

For a no-evolution model for the galaxy population, the 



Figure 9. The Bj band correlation function amplitude. As with 
the Rf band data, the amplitude of the correlation function is 
significantly weaker than the / band correlation function. At mag- 
nitudes fainter than Bj ^ 22, the correlation function rapidly 
decreases with limiting magnitude. There is also a decrease in 
the discrepancy between the SGP and F855 measurements with 
increasing limiting magnitude. 
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Figure 10. The U band correlation function amplitude. The am- 
plitude of the correlation function is consistently weaker than the 
fit to the I band data shown with the dashed line. 



blue bands will sample galaxies with bluer colours due to 
the large fc-corrections of early-type galaxies (Coleman, Wu 
& Weedman 1980) . As shown in Figure |ll|, the local popula- 
tion of blue {U — Bj ^ 0) galaxies is dominated by late-type 
galaxies. Measurements of the galaxy clustering in the lo- 
cal universe with morphology-selected catalogues show the 
clustering of late-type galaxies is considerably weaker than 
early-type galaxies (Davis & Geller 1976, Loveday et al. 
1995). If the colours and clustering of 2 ~ 0.4 late and early 
type galaxies are similar to 2 ~ galaxies, the decrease in 
the amplitude of the Bj angular correlation function could 
be a selection effect. While it is impossible to determine the 
morphologies of the Bj > 21 galaxies with this catalogue, 
it should be possible to use colour selection to select early 
and late type galaxies over a range of redshifts. 
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Figure 11. The colours of RC3 catalogue (de Vaucouleurs et 
al. 1991) galaxies by morphological type. Photometry is from 
Prugniel & Heraudeau (1998) while morphological classification 
of galaxies is from de Vaucouleurs et al. (1991). The correlation 
between morphology and colour can be clearly seen with late type 
galaxies generally having bluer colours than early types. 

7 THE CORRELATION FUNCTION OF 
COLOUR SELECTED GALAXIES 

The obvious colour selection criteria for galaxies is a single 
colour cut in the deepest bands available {Bj and Rf). How- 
ever, as shown in Figure the colour of individual galaxy 
types varies with redshift. Even with a single colour cut, 
the blue subsample will generally select later type galaxies 
than the red subsample and blue subsamples generally show 
weak clustering (Infante & Pritchet 1995, Roche et al. 1996). 
A Bj — Rf cut that varies with magnitude may select the 
same population over a range of redshifts but the selection 
criteria would depend on the cosmological parameters used. 
Also, colour selection of a fraction of galaxies (i.e. the red- 
dest 20% of the catalogue) may not be effective due to the 
changing morphological mix of galaxies with magnitude. 

Selection of galaxies with two or more colour criteria 
should allow the selection of galaxy types over a large range 
of redshifts without strong dependence on cosmological pa- 
rameters. Figure |l2| shows the predicted colours of early, 
late and irregular type < 2; < 0.8 galaxies from Fukugita, 
Shimasaku & Ichikawa (1995). Colour selection criteria for 
a red and blue subsample are also shown. The samples are 
limited to Bj < 21.6 due to the U — 22 magnitude limit 
of the catalogues and the U — B,j > 0.4 selection criteria. 
As most Bj < 21.6 galaxies are at redshifts z < 0.3, most 
galaxies in the subsamples are selected with the U — Bj 
selection criteria. Galaxies with U — Bj lower limits and 
Bj — Rf upper limits have been included in the subsamples 
to prevent incompleteness. Figure ^ shows the colours of 
galaxies in the SGP with the subsample selection criteria. 
The location of the Bj < 19.5 galaxy locus is similar to that 
for RC3 catalogue galaxies in Figure |ll| though at Bj > 19.5 
there is an increasing fraction of very blue galaxies. Galaxy 
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Figure 12. The colour selection criteria for the red and blue 
subsamples. Colours of early, late and irregular type galaxies at 
z = 0, z = 0.2, z = 0.5 and z = 0.8 from Fukugita, Shimasaku & 
Ichikawa (1995) are shown. The colour selection criteria for each 
subsample are shown with the dot-dash line. The blue subsample 
consists of galaxies below and to the left of the line while the 
red subsample consists of galaxies above and to the right of the 
line. The blue subsample is dominated by late and irregular type 
galaxies while the red subsample selects early type galaxies. 



number counts for the 2 subsamples are shown in Figure |l4|. 
While the SGP does include significantly more clusters than 
F855 at low redshift, there is good agreement between the 
number counts across the magnitude range observed. 

While the selection criteria are relatively simple, it can 
be clearly seen that the red and blue subsample should select 
early and late type galaxies respectively. The blue subsample 
should also contain most of the galaxies with significant star 
formation rates while the red subsample should contain more 
passive galaxies. As shown in Figure ^ comparison with the 
spectroscopic catalogue of CoUess et al. (1990) shows that 
most of the [Oil] 3727A emitters detected in the SGP and 
F855 are included in the blue subsample. To measure the 
clustering of starforming galaxies Cole et al. (1994) selected 
galaxies with [Oil] 3727A equivalent widths greater than 
I9A; all galaxies matching this selection criteria would be 
included in the blue subsample for both fields. Interestingly, 
the galaxies without [Oil] 3727A emission do not show such 
an obvious trend with similar numbers in both subsamples. 
However, the small number of galaxies without [OH] 3727A 
emission results in these galaxies comprising less than 15% 
of the total of blue galaxies. 

The Bj < 21.5 angular correlation functions of the blue 
and red subsamples of F855 are shown in Figure ^ Vi- 
sual inspection shows the significant difference in clustering 
strength and the value of 7 for the two subsamples with the 
red subsample being strongly clustered and having a higher 
value of 7. This trend is consistent with measurements of the 
angular correlation functions of morphology selected cata- 
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Figure 13. Galaxy colour-colour diagrams for the SGP shown with the selection criteria for the two subsamples. 
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Figure 14. A plot of the red and blue subsample number counts 
for the SGP and F855. 



logues where 7 ~ 1.8 for early-type galaxies and 7 ~ 1.5 for 
late- type galaxies (Loveday et al. 1995). 

Estimates of 7 and as a function of limiting mag- 

nitude are listed in Table W. The z ~ 0.11 clusters in the 
SGP have been retained in the data as there is no obvious 
justification for rejecting them from the sample. While 7 dif- 
fers significantly between the two subsamples, the values of 
7 are within Icr of being constant as a function of limiting 
magnitude. The values of uj(1') for the blue subsamples are 
remarkably similar for the SGP and F855. However, the red 
subsamples show large differences with the amplitude of the 
clustering in the SGP being ~ 100% larger than F855. 

To model the spatial correlation function, a model of 



m 

I o 

ZD 



1 X 


1 1 1 1 1 1 1 1 1 1 
Red Subsample 


- n 


t . no [Oil] - 
t X [Oil] <19A - 
[Oil] >19A ; 


X 
X 



HP 




: » ° . 




X 

• 

" Blue Subsample 
, , , , 


1 , , , , 1 , , , , 



CM 
I 



Bj-Rf 

Figure 15. A colour-colour diagram of Bj < 21.5 galaxies in the 
SGP and F855 with spectroscopic observations by Colless et al. 
(1990). Lower limits va U — Bj are shown with arrows. A trend 
towards bluer colours with [Oil] 3727A emission can be observed 
with strong [Oil] 3727A emitters restricted to the blue subsample. 



the redshift distribution for the red and blue subsamples is 
required. Luminosity functions for galaxies selected by mor- 
phology are available but it is not clear if luminosity is more 
or less strongly correlated with colour than morphology. We 
therefore use the U — Bj < 0.2 and U — Bj > 0.2 galaxy 
luminosity functions of Metcalfe et al. (1998). Approximate 
fc— corrections of k{z) = 3z and k{z) = z are used for the red 
and blue subsample respectively. We use the Metcalfe et al. 
(1998) luminosity functions rather than the redshift distri- 
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Table 4. Measured parameters for t<j{9) for the red and blue subsamples determined from power law fits to data between 10" and 0.3°. 
Field SGP F855 

Sample Magnitude Range Ng^i 7 <.^(1') X 10^ Ngai 7 ^^(1') X 10^ 

Blue 18.0 < Bj < 20.0 2963 1.45 ±0.51 212 ± 3200 3354 1.41 ± 0.58 170 ± 2710 

Blue 18.0 < Bj < 21.0 9079 1.49 ±0.18 138 ± 27 9228 1.61 ±0.38 110 ±39 

Blue 18.0 < Bj < 21.5 15579 1.57 ±0.20 86 ± 19 14329 1.51 ±0.32 75 ± 32 

Red 18.0 < Bj < 20.0 1993 1.65 ±0.22 1339 ± 278 1851 2.09 ±0.17 954 ± 210 

Red 18.0 < Bj < 21.0 5552 1.69 ±0.12 782 ± 73 5424 1.89 ±0.13 480 ± 82 

Red 18.0 < < 21.5 9317 1.72 ±0.08 616 ± 51 9922 1.83 ±0.14 319 ± 44 



Bj<21.5 Galaxy Angular Correlation Function 




6 (Degrees) 



Figure 16. The Bj < 21.5 angular correlation function for 
the red and blue subsamples of F855. The red subsample shows 
stronger clustering and a higher value of 7, consistent with the 
strong clustering of elliptical and early type galaxies in the local 
universe. 



bution model as the steep slope of faint end of the luminosity 
functions results in skewed redshift distributions. 

Plots of the estimated median redshift for the red and 
blue luminosity functions are shown in Figures |l^ and |l8[ 
The SGP and/or F855 overlap the Bj magnitude limited 
redshift surveys of CoUess et al. (1990) and Ratcliffe et al. 
(1998) which have been used to measure the median red- 
shifts of the red and blue subsamples as a function of mag- 
nitude. In addition, ~ 700 galaxy redshifts are available from 
the NASA/IPAC Extragalactic Database (NED) and these 
have been used to show the approximate median redshift as 
a function of magnitude. For both subsamples, the model is 
a reasonable estimate of the redshift with few data points 
being more than 20% from the model redshift estimate. 

Figure |l^ shows the amplitude of the SGP and F855 
red subsamples as a function of limiting magnitude with 7 
fixed at 1.8. A model of the early type galaxy correlation 
function with ro = 5.9/i~^Mpc (Loveday et al. 1995) and 
clustering fixed in physical coordinates (e — 0) is also shown. 
Clustering fixed in physical coordinates would be applicable 
to galaxies in gravitationally bound clusters. The model, 
shown with the dashed line, does not fit either set of data but 
this is not unexpected as the SGP and F855 fields exhibit 
significantly different clustering amplitudes. However, the 
clustering in both fields is consistently stronger than the 



Red Subsample 

in 




o I I < I 1 I 1 I 

o 16 18 20 22 



Bj Magnitude 

Figure 17. The median redshift as a function of magnitude for 
the red subsample. Data points from Bj band magnitude limited 
surveys are shown with circles while data using NED redshifts is 
shown with dots. The median redshift estimate and data points 
have been determined using 0.5 magnitude wide bins. As the es- 
timates of the correlation function use wider magnitude bins, the 
estimates of the median redshift in this plot are higher than in 
Figure ha. 



clustering of all Bj galaxies in F855. This is consistent with 
the red subsample being dominated by strongly clustered 
early type galaxies. 

Figure bd plots the amplitude of the blue subsample cor- 
relation function derived from power law fits with 7 fixed 
to 1.5. Unlike the red subsample, the amplitudes derived 
from the SGP and F855 data agree at all magnitudes. If 
dwarf galaxies are increasing the amplitude of the SGP 
correlation function, this implies that they are red galax- 
ies such as dwarf ellipticals within the clusters. A model 
with extremely weak clustering (ro — 2/!,~^Mpc) fixed in co- 
moving coordinates is also shown in Figure The ampli- 
tude of the clustering is similar to that estimated for faint 
galaxies by Efstathiou (1995) and z ~ 0.5 galaxies by Le 
Fevre et al. (1996). The model slightly underestimates the 
strength of the clustering and a model with stronger clus- 
tering, ro = 3/i~^Mpc, is a better fit to the data. This is 
significantly weaker than ro ~ 4.4/i~^Mpc, the measured 
clustering of late type galaxies in the local Universe (Love- 
day et al. 1995). The only large sample of low z galaxies with 
similar properties (ro ~ 3^~^Mpc, 7 ~ 1.8) are galaxies 
with [Oil] or Ha emission lines with large equivalent widths 
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Figure 18. The median redsiiift as a function of magnitude for 
the blue subsample. Data points from Bj band magnitude limited 
surveys are shown with circles while data using NED redshifts is 
shown with dots is shown with dots. The median redshift esti- 
mate and data points have been determined using 0.5 magnitude 
wide bins. As the estimates of the correlation function use wider 
magnitude bins, the estimates of the median redshift in this plot 
are higher than in Figure 201 



Blue Subsample — 
n=0.2, A=0.0, Ho=75.0 kms"'Mpc"' 
Model r|j=2.0 h"' Mpc, 7=1.50. e=-1.50 
Model r„=3.0 h"' Mpc, 7=1.50. e=-1.50 
Power Law Fit to SGP Data with Ict Error Bars 
Power Low Fit to F855 Data with Iff Error Bars 
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Figure 20. The amplitude of the clustering for the blue sub- 
samples as a function of magnitude. The value of 7 has been 
fixed to 1.5. The blue galaxies are extremely weakly clustered 
with a similar spatial correlation function to faint Bj galaxies. 
As the Bj ~ 23 number counts are dominated by blue galaxies 
(Guhathakurta, Tyson & Majewski 1990), this is strong evidence 
for star forming galaxies being weakly clustered from 2 ~ 0.4 
until the present epoch. 
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Figure 19. The amplitude of the clustering of the red subsamples 
as a function of magnitude. The value of 7 has been fixed to 1.8. 
The clustering strength in the two fields differs by ~ 100% at all 
magnitudes. Both fields show strong clustering, consistent with 
the red subsamples being dominated by early type galaxies. 



from the Stromlo-APM survey (Loveday, Tresse & Maddox 
1999). However, the Stromlo-APM sample has significantly 
stronger clustering than the blue subsample or faint blue 
galaxies on scales ^ 2/i~^Mpc. 

An underestimate of the redshift by <: 50% could ex- 
plain the weak clustering at Bj ~ 20 but this would be diffi- 
cult to reconcile with the redshift data in Figure Altering 
the assumed cosmological parameters can change estimates 
of spatial correlation function but the effect is negligible at 
~ and is less than 30% ai z — 0.2. A plausible ex- 
planation is that the clustering of galaxies is more strongly 
correlated with colour and stellar population than morphol- 
ogy. This would explain why no local population of galaxies 
selected by morphology displays the weak clustering of faint 



blue galaxies. This may also be the first detection of large 
population of galaxies at low redshift with similar clustering 
properties to faint blue galaxies. 



8 SUMMARY 

We have used deep multicolour galaxy catalogues of two 
5° X 5° fields to study clustering from 2 ~ to 2: ~ 0.4. The 
key conclusions of this paper are: 

(i) The galaxy spatial correlation function is a power law 
on comoving scales less than 15/i~^Mpc. At larger scales, the 
correlation function is consistent with a power law though 
a break in the correlation function is not inconsistent with 
the data. 

(ii) Despite the large fields-of-view, there are significant 
differences in the measured amplitude of the clustering; with 
the possible exception of blue galaxies. It is clear that fields 
larger than 100n° are required to accurately measure the 
clustering of Bj ~ 22 galaxies. 

(iii) Dwarf galaxies in relatively nearby clusters (2 ~ 
0.11) may effect estimates of faint galaxy correlation func- 
tion. The effect is colour dependent with the clustering of red 
galaxies varying significantly between the 2 fields observed. 

(iv) The clustering properties of galaxies strongly de- 
pend on the band used to select the catalogue. Bluer bands 
show weaker clustering than red bands and there is a rapid 
decline of the amplitude of the Bj correlation function at 
faint magnitudes. It is probably inappropriate to fit a simple 
clustering model to correlation functions derived from single 
band imaging due to the changing morphological mix with 
magnitude. 

(v) The clustering properties of galaxies strongly de- 
pend on colour. Such behaviour is consistent with colour be- 
ing correlated with morphological type. Red galaxies (early 
types) exhibit stronger clustering with larger values of 7 
than blue galaxies (late and irregular types). 
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(vi) Blue galaxies have extremely weak clustering with 
ro ^ 3^~^Mpc. This is considerably weaker than the cluster- 
ing of late type galaxies and is consistent with the clustering 
of galaxies being more strongly correlated with colour and 
stellar population than morphology. 

(vii) The clustering of Bj < 21.5 blue galaxies is com- 
parable to Bj > 23 blue galaxies. This is strong evidence for 
star forming galaxies being weakly clustered from z ~ 0.4 
until the present epoch. 
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